OSCILLATION OF LINEAR ORDINARY DIFFERENTIAL 

EQUATIONS: 
ON A THEOREM BY A. GRIGORIEV 

SERGEI YAKOVENKO* 



Abstract. We give a simplified proof and an improvement of a recent 
theorem by A. Grigoriev, placing an upper bound for the number of 
roots of linear combinations of solutions to systems of linear equations 
with polynomial or rational coefficients. 



1. Background on counting zeros of solutions of linear 
ordinary differential equations 

1.1. De la Valee Poussin theorem and Novikov's counterexample. 

A linear nth order homogeneous differential equation 

+ ai(t) + • • • + a„_i(t) y' + a„(t) y = 0, y^'^) = (1.1) 

with real continuous coefficients aj{t) is called disconjugate {Chebyshev, non- 
oscillating) on a real interval [a, f3] C M, if any nontrivial solution y{t) of 
this equation has at most n — 1 isolated root on this interval. 

A theorem by C. de la Vallee Poussin |dlVP29j asserts that any equation 
(|l.lj) is disconjugate on any interval sufficiently short relative to the mag- 
nitude of the coefficients of the equation. More precisely, if |aj(t)| ^ bj on 
[a, f5] and Yl^=i ~ o^P/j- < 1) then 1)1.1(1 is disconjugate. This allows to 
place a rather accurate upper bound on the number of isolated roots of any 
solution of a known differential equation in terms of the length of the inter- 
val on which the solution is considered and the magnitude of the coefficients 
of the equation. A complex analog of this theorem was obtained in |Yak99j 
for linear homogeneous nth order equations with holomorphic coefficients in 
a polygonal complex domain t S [/ d C. 
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x{t) = A{t)x{t), 
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A natural generalization of these results for systems of first order linear 
homogeneous differential equations of the form 

•. ; , (1-2) 

with variable coefficients matrix A{t), might deal with a uniform upper 
bound for the number of isolated zeros of an arbitrary linear combination 
of components y{t) = ci xi{t) + • • • + CnXn{t), in terms of the length of 
the interval t G [a, f3] where the solution is considered, and the magnitude 
max^gjQ,^^] of the matrix norm on this interval. Among other reasons, 

such result could be expected since for any such linear combination one can 
explicitly derive a linear nth order equation of the form (jLlj) , as explained 

in 32 

Unfortunately, nothing like this is possible. In |Nov01j D. Novikov con- 
structed a very simple system with the following properties: 

(1) the coefficient matrix A{t) is a real matrix polynomial of degree m, 

m—l 

A{t) = Akt'', Afc G Mat„(C), k = 0,...,m-l, (1.3) 

j=0 

in the minimal dimension n = 2, 

(2) the matrix norm ^ 1 everywhere in the disk < 2} C C; 

(3) the first component xi{t) of some solution x{t) of the respective 
system (jl.2|) has as many isolated roots in the unit disk {\t\ < 1} as 
necessary, provided that the degree m of the matrix polynomial can 
be arbitrarily large. 

When reducing the Novikov's system to a second order equation (jLlj) . 
one encounters uncontrollably large rational coefficients aj{t) = Pj{t) /po{t): 
while the numerators Pj{t) can be explicitly majorized, the denominator 
Po(i) can be arbitrarily small (either coefficients- wise, or uniformly on large 
intervals) . 

1.2. Bounded meandering principle. In an attempt to circumvent the 
appearance of the small denominators, one can declare all entries of the ma- 
trix coefficients ^i) • • • i Am-i of the polynomial ()1.3() to be new phase 

2 

variables A G M" ™ playing the same role as the initial variables x £ W^. 
Then the system ()1.3() coupled with the trivial equations A = 0, t = 1, 
becomes nonlinear, but all numeric coefficients of this nonlinear system be- 
come integer (in fact, only zeros or ones). This lattice structure allows to 
avoid appearance of the small denominators (any integer number is either 
zero or at least one in the absolute value), but the "reduction" to a scalar 
linear equation becomes considerably more complicated. 

In |JN Y99j the problem of counting roots of solutions of systems of poly- 
nomial differential equations was reduced to estimating the length of an 
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ascending chain of polynomial ideals. The ultimate result established in 
|NY99j is a general "bounded meandering principle": for any system of 
polynomial ordinary differential equations and any polynomial combination 
of components of its solutions, the number of isolated roots of this combi- 
nation on any finite interval can be explicitly bounded in terms the degrees 
of all polynomials, dimension of the system, the length of the interval and 
the magnitude of the coefficients of the equations. More specifically, if the 
coefficients of the system of n differential equations and the length of the 
interval on which the combination is considered, do not exceed some upper 
bound M > 2, and the degrees of all polynomials do not exceed m, then the 
number of isolated roots of any polynomial combination of solutions does 
not exceed M'^, where q = q{n, m) is an explicit computable function of two 
natural arguments, growing no faster than the iterated exponent of order 4: 



where 0(1) is an absolute constant |NY99j . 

1.3. Grigoriev theorem. The approach based on the bounded meander- 
ing principle, completely ignores the linearity of the initial system (|1.3j) : 
the bound for the number of isolated zeros is given by the same extremely 
excessive expression (|1.4j) . 

In 2001 in his Ph. D. thesis [(xriOlj recently published in the revised form 
[(jriOSj , Alexei Grigoriev returned to the study of the system (|1.3j) by "purely 
linear" methods. The principal obstruction on this way, already mentioned 
above, is occurrence of the "singular degeneracy" : for certain values of the 
parameters A, the resulting equation has the leading coefficient (before the 
highest derivative) vanishing identically. In the classical language, the para- 
metric family of differential equations near such values of the parameters 
would exhibit a singular perturbation. 

Appearance of singular perturbations can indeed be a source of highly 
oscillating solutions. The simplest examples can be easily constructed al- 
ready in the class of linear equations with constant coefficients. Yet in the 
situation when the equation is obtained from a system ()1.3|) . we have an 
additional bit of information: all solutions of the derived equation depend 
analytically on the parameters. 

In jYak99l p. 520] it was conjectured that even in the family exhibiting 
a singular perturbation, all solutions which depend analytically on the pa- 
rameters admit an effective uniform upper bound for the number of isolated 
zeros. A. Grigoriev proved a similar conjecture for the equations with ra- 
tional coefficients whose all solutions remain analytic, and achieved a much 
better upper bound than the tower ()1.4|) . In the simplest settings his re- 
sult can be formulated for a system of the form H1.3() with the polynomial 
coefficients. 

Being polynomial, such system can be considered also for complex values 
of time t. We are interested in an upper bound for the number of isolated 




^24n lnm + 0(l) 



(1.4) 
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roots of any linear combination in a bounded domain. Without loss of 
generality we may assume (rescaling the variable t, if necessary) that this 
domain is the unit disk D = < 1}, while all matrix coefficients are 
bounded by the constant M: 

\t\<l, \Ak\i^M, k = 0,...,m-l. (1.5) 

The following result can be obtained by slightly improving the arguments 
by A. Grigoriev. 

Theorem 1. The number of isolated zeros of any linear combination of 
components of solutions of the system constrained by the restrictions 

H1.5() in the unit disk does not exceed 

[max(M,2)]2 (1.6) 
where the constant in O is absolute. 

Remark 1. The bound H1.6() seems to be still rather excessive. For instance, 
using the standard arguments relating the growth of analytic functions with 
the number of their zeros via the Jensen inequality |Lev80[ IIY96j , we can 
prove that for most linear combinations, the number of zeros does not exceed 
we consider fixed parameter. This is a far shot from the 

double exponential bound (|1.6|1 . 

Remark 2. Grigoriev in [(irrin.Sj formulates this theorem only for the case 
M = 1, and only as a remark; besides, the bound that he claims in this 

case, is somewhat worse: 22'"™' instead of 22"'" ™' in (jOni . 

The goal of this work is to reveal some geometric ideas lying behind 
the Grigoriev's proof which (despite the recent revision) remains still highly 
technical and sometimes rather obscure. Besides, we discuss how some other 
classes of differential equations can be treated in the similar way, in partic- 
ular, why the hypergeometric systems are better than a general Fuchsian 
system from the point of view of counting the roots. This is especially im- 
portant in connection with the tangential Hilbert problem on the number 
of zeros of Abelian integrals |NYni[ lYakOlj . 

2. From system to a high order equation 

2.1. Variations on the theme of de la Vallee Poussin. In what follows 
we will use the £^-norm on the space of polynomials in one or several vari- 
ables: by definition, the norm of a polynomial is the sum of absolute values 
of all (real or complex) coefficients, 

/gC[x] =C[xi,...,:c„], /=J^c„x", ll/llcN = J^lcal- (2.1) 

The main advantage of this norm is its multiplicativity: ^ ||/|| 

Consider the equation (ll.lj) with polynomial coefficients ao,ai, . . . ,a„ G 
C[t] of known degrees. Note that, since the equation is homogeneous, the 
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coefficients are defined modulo a common factor. The following result is a 
variation on the theme of Theorem 2' from IY96 : if the principal coeffi- 
cient before the highest derivative is not too uniformly small while all other 
coefficients explicitly bounded from above, then the number of zeros of any 
holomorphic solution is explicitly bounded (non-holomorphic solutions with 
singularities may disobey this principle). 

Lemma 1. Assume that the coefficients of the equation of order n 

are polynomials of degree ^ d, and their norms \\aj\\ = ||oj||c[t] satisfy the 
restrictions 

||ao|| = 1, llofcll ^ A; = l,...,n. (2-2) 

Then any entire solution to (|1.H) may have at most Kn ■ 2^^'^^ isolated 
roots in the unit disk D = < 1} C C, where the constant in 0(d) is 
absolute. 

Proof. We start with observing that a polynomial gq G C[t] of degree ^ d 
normalized by the condition ||oo|| = 1, cannot be uniformly small on, say, 
the annulus {1 < |t| < 2}. More precisely, there exists a circle {\t\ = r}, 
1 < r < 2, on which |ao(*)| is bounded from below, |ao(t)| ^ 2~^^'^\ 

Indeed, the condition ||ao|| = 1 means that one of the coefficients of the 
polynomial ao is at least l/d in the absolute value. Applying the Cauchy 
formula for the respective derivative of ag, we conclude that 



1/d ^ 



2ni s^'+i 



^ max |ao(s) 

|s|=l 



Therefore that there exists a point sq on the unit circle \t\ = 1, on which 
|ao(so)| is at least 1/d. Applying |Lev8n[ Theorem 4, p. 79] to the polynomial 
(^ois) = oo(s — so)/d, we conclude that the lower bound |ao(s)| ^ 2~'^('^) with 
some absolute constant in 0{d) holds on the complement to a finite union 
of disks of the total diameter less than 1. Such disks cannot intersect all 
concentric circles {\t\ = r} , 1 < r < 2, which proves the required assertion. 

On this circle {\t\ = r} all other coefficients are bounded from above, 
|aj(t)| ^ K ■ 2*^('^) for all j = 1, . . . , n. Restricting the equation on the 

circle parameterized by t = rexpis, s G [0, 27r], and applying Corollary 2.7 
from |Yak99[ p. 507], we conclude that the variation of argument of any 
solution y{t) of the equation on this circular segment, can be at most 

Kn-2C'W. 

If y{t) is holomorphic in the disk {\t\ < r}, in particular, if this solution 
is an entire function, then by the argument principle the number of its zeros 
does in the smaller disk D does not exceed Kn ■ 2*^^*^^ as asserted. □ 

Remark 3. The requirement that the solution is entire, is excessive. In fact, 
it is sufficient that the solution of the equation is analytic in a disk of any 
finite radius larger than 1. For instance, the proof below works literally if 
y(-) is analytic in the disk of radius 6. If necessary, the number 6 can be 
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replaced with only minor modifications by any number bigger than 1, but 
this choice will affect the constant in 0{d). 

2.2. Derivation of the principal equation. Consider now the linear sys- 
tem We will denote by A = {Xijk} G C" the collection of all entries 
of all matrices A^, so that the coefficients matrix of the system would 
be a "universal" matrix polynomial A{t,\) from the ring A], whose en- 
tries have norm m in this ring. 

Without loss of generality, when proving Theorem Q it is sufficient to 
consider only the first coordinate y{t) = xi{t) of an arbitrary solution of 

jni). 

The formulas for higher derivatives of the function y(-) can also be ex- 
pressed by linear combination of the components A), . . . , Xn{t, A) with 
polynomial coefficients. Indeed, assume inductively that 

2/^'^(t) = efc(i,A)-x(t,A), 

where A), = 0, 1, . . . , n is a polynomial (row) vector function. Then 
after derivation by virtue of the system (|1.2() we obtain 

y{k+i) {t) = ^. x{t, A) + 6 • A{t, A) • x{t, A) = + CkA^ ■ x{t, A). 

This allows to define the vector functions starting from from the initial 
vector-function ^o(i, A) = (1, 0, . . . , 0) G C". 

Consider the collection G C[t, A] (8) C" of polynomial vector-functions 
defined by the recursive formula 

^k+i = ^^^k+^k-A{t,X), k = 0,l,...,n-l. (2.3) 

The process (|2.3p is so explicit that the following estimates are almost 
obvious. 

Lemma 2. The collection of polynomial vector functions S,o,---,£,n G 
Z[t,A] C" satisfies the following inequalities: 

(1) the vector polynomials (,o, ■ ■ ■ cli"^ integral, i.e., they belong to the 
lattice A] (g)z c C[t, A] ®c C", 

(2) deg.^fc ^ km for all k = 0, . . . ,n, 

(3) the wedge product of any n vectors out of this collection is a scalar 
integral polynomial from Z[t,X] of the degree at most ^n{n — l)m 
and the norm not exceeding n\ (2nm)" , 

(4) the wedge product A • • • A G Z[t, A] is not identically zero in 
Z[t,A]. 

Proof. Everywhere until the end of the proof we abbreviate deg = degj )^ 
and II • II = II • ||c[t,A]- 

The fact that all coefficients of are integral over t. A, follows immedi- 
ately from the fact that the "universal" matrix polynomial A{t, A) is integral, 
A G Z[t,A] ®z Mat(n,Z). 
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The inequality 

deg Cfc+i ^ deg + deg A{t, A) ^ deg + m 

proves the assertion about the degrees of the functions themselves. The 
Laplace formula for the complete expansion of the (nxn)-determinant proves 
the assertion for the degrees of the wedge products. 

The assertion about the norm follows from the recursive inequality 

Uk+i\\ ^ (degCfc + ^ 2nm||^fc||, 

where by the norm of the polynomial i^ector-function we understand the 
maximum of the norms of its scalar polynomial components. Therefore 
each entry of the (n x n)-determinant has the norm not exceeding (2nm)", 
and by the Laplace expansion we have the sum of n! terms, each of the norm 
at most [(2nm)"]" by the multiplicativity of the norm. 

To prove the last assertion on the nondegeneracy, it is sufficient to note 
that for the specific value of the variables A corresponding to the con- 
stant" (i.e., independent of t) polynomial A = Aq without higher terms, 
the vectors = ^o^o are also "constant" vectors. For the generic choice 
of Aq E Mat(n, C) they are obviously linear independent over C for all 
k = 0,1, . . . ,n — 1, and hence the wedge product is not identically zero. □ 

Corollary 3. The first coordinate of any solution of the system is an 

entire function satisfying a linear nth order differential equation of the form 

ao{t. A) y(") + ai{t. A) y^""!) + • • • + a„(t, A) y = 0, (2.4) 
for the corresponding value of the parameters A. 

The coefficients aj{t,\) € Z[t,A] of this equation satisfy the following 
conditions: 

(1) all aj £ 7j[t,X] are integral (lattice) polynomials, 

(2) degixaj ^ ^ri{n - l)m, 

(3) ||«j||c[t,A] ^ nl (2nm)"^ 

(4) the leading coefficient ao{t,X) does not vanish identically. 

Proof. For arguments of pure dimensionality, the vectors Co, £ 
Z[t, A] Z" C C(t, A) C must be linear dependent over the field of ra- 
tional functions C(t,A), 

aoCn H h ttnCo = 0. 

This identity implies a linear equation of order n for the function y{t) = C,o-x. 
The coefficients of this dependence can be found (by the Cramer rule) as 
the ratios of the wedge products of the vectors of degree n. These were 
already described in Lemma El D 
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2.3. Investigation of the principal equation. The fact that the leading 
coefficient ao of the equation (|2.4|) does not vanish identically in A], does 
not imply that for some values of the parameters A it becomes a zero polyno- 
mial of t. Such values of the parameters correspond to degenerated systems 
for which the linear dependence over the field C(t) between derivatives of the 
first coordinate occurs earlier than on the nth step. In classical terminology, 
this phenomenon should be described as a singular perturbation, when the 
leading coefficient of a differential equation vanishes identically. 

The main result by Grigoriev, which was conjectured in |Yak99| §4.5, 
p. 520], is the assertion that in fact this singular perturbation is only ap- 
parent, namely, that for these exceptional values of the parameter, all other 
coefficients of the principal equation 1)2. 4|) are also vanishing identically, and 
their ratios are explicitly bounded for all nearby values of the parameters A. 

More precisely, consider the functions 

6,- : C^'"* ^ M+, b,{X) = \\a,{;X)\\c[t], i = 0,l,...,n, (2.5) 

the norms of the coefficients aj{t,X) of the equation ()2.4() considered as 
univariate polynomials. The zero locus U of the function bo corresponds to 
the "singular degeneracy" of this equation. 

The cornerstone of the proof of Theorem^ is the following analytic result 
discussed in details and proved in the next section ^ 

Lemma 4 (Principal Lemma). The ratios bj/b^ are locally bounded on 
C^^""^!! for all j = 1,... ,n. 

In the remaining part of this section we derive Theorem ^ from Lemma HI 
The first step is the following Corollary, transforming the quantitative state- 
ment of Lemma m into a qualitative bound. 

Lemma 5. In the notations (|2.4j) - (|2.5j) . 

^ 2)'"'"''"' ' J = 1, • • • , n, (2.6) 

where the supremum is taken over all values X ^ U = {60(A) = 0}. 

Proof of the Lemma. Consider instead of the norms bj{X) = ||oj(-, A)||c[f] 
equal to the sum of modules coefficients of the univariate polynomials 
aj(-,A), the functions Cj(A) equal to the sum of the squares of absolute 
of those coefficients. The advantage is that the functions Cj{X) are polyno- 
mials in Re A and ImA. On the other hand, by the Cauchy inequality we 
have two-sided equivalence, 

Cj{X) ^ b]{X) ^ ncj{X), 

therefore it will be sufficient to prove (|2.6j) for the ratios Cj{X)/co{X): the 
additional power 1/2 and the factor y/n will be obviously absorbed in the 
O-symbol. 
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Without loss of generality we may assume that the number M is natural. 
In the real space M^" x M+ (the first factor is the realification of C" 
consider the subset 

Sj,M = {{X,u): co(A) ^0, u^O, |A| ^ M, co(A)m ^ Cj{X)}. 

This subset is semialgebraic (defined by polynomial equalities and inequali- 
ties) . The coefficients of these equalities and inequalities are integer numbers 
which do not exceed max.{M,N), where N = re!(2nm)"' , and the degrees 
do not exceed v?m. 

In addition, this set is bounded by Lemma 01 Sj^m is a subgraph of a 
locally bounded function defined on a bounded set {|A| ^ M}. 

A bounded semialgebraic set defined by integer (lattice) polynomial 
(in)equalities of degree 6 and the norms at most /x in M^, admits an 
explicit, effective and rather accurate upper estimate of its size: by 
|(j;V881 IHHSQOl |BPR03j, such set always belongs to the cube of the size 
l/^^"^ centered at the origin in M*^. 

Substituting the known values of the norms and the degrees, we obtain 
the required upper bound for the u-coordinate, as required. Note the most 
crucial parameter is the dimension v equal in our case to n^m (the total 
number of different parameters describing the "universal" system 1)1.3^ ) . All 
other dependencies on n, m are absorbed into the term O(n^m) occurring 
at the second floor of the tower function. □ 

Example 1. Consider the real rational function of two variables r{x,y) = 
(ax^ + by'^)/ {cx^ + by^). If the coefficients a, 6, c, d G M are all positive and 
ad — be ^ 0, then this function is locally bounded near (x,y) = (0,0). Yet 
the limits along different rays y = ux, u G M, are different and vary between 
a/c and b/d. Knowing only the upper bound // for the coefficients a,b,c,d 
yields no local upper bound for r{x, y), but if these coefficients are integral, 
then fj,~^ ^ r{x,y) ^ fi. 

2.4. Proof of Theorem ^ For all values of the parameters A outside the 
degeneracy locus U, the assertion of Theorem ^ follows from Lemma ^ into 

which one should substitute the value of K = max(M, 2)^*^'" given by the 
LemmaEl For X £ U the principal equation (|2.4|) degenerates into the trivial 
identity = 0, but solutions (the first coordinates of solutions of the system 
(|l..Sjl ) remain entire functions of t and A. The number of isolated zeros of 
an analytic function in the open disk < 1} depends semi-continuously 
on the parameters and hence the expression given in the Theorem, serves 
as an upper bound for the number of isolated zeros also for the exceptional 
values of the parameters A E 17. □ 

3. Analytic one-parameter families of linear subspaces 

The proof of the Principal Lemma^Jis based on the following observation: 
a linear space spanned by any number of vectors analytically depending on 
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one parameter, depends very regularly on this parameter even if the vectors 
exhibit extra linear dependence for isolated values of the parameter. 

3.1. Finite-dimensional case. We first illustrate the above observation 
for subspaces of a finite-dimensional space. Everywhere below the parameter 
e is one-dimensional (real or complex) and local, varying in a neighborhood 
of zero. 

Proposition 6. Let vi, . . . ,Vm- (M"*^, 0) M" be any real analytic germs of 
vector-functions, m ^ n. Then there exist k real analytic germs of vector 
functions wi, . . . ,Wk- (M-^, 0) — > M", k ^ n, which have the constant rank k 
and span the same linear subspace in M" as the vectors fi(e), . . . ,fm(£) for 
all e^O. 

Proof. The proof goes by induction in m. If m = 1 and vi{e) = 0, then 
there is nothing to prove, otherwise ui(e) = e'^wi{e) with wi an analytic 
vector-function such that wi{0) ^ 0. Clearly, the linear spaces spanned by 
vi{e) and wi(e) coincide for e 7^ 0. 

To prove the induction step, assume that the vectors wi{e), . . . ,Wk-i{£) 
are linear independent and span the same linear space L{e) as the vectors 
vi{e), . . . ,Vm-ii£) for all e 7^ 0. If Vmi^) £ L{e) for all e, then vji, . . . ,Wk-i 
span the same space as vi, . . . ,Vm- Otherwise choose any linear subspace K 
in complementary to L{0) (hence to all L(e) for small e) and the projec- 
tion TTg : — > L{e) parallel to K. The vector function fm(e) — TTsVmis)- By 
assumption, it is not identically zero and is real analytic, since vr^ depends 
analytically on e by the Implicit function theorem. Therefore, 

for some finite power v. The vector Wk{0) is transversal to L{0), hence the 
system wi{e),. . . ,Wk{s) is linear independent for all small e and generate 
the same space as vi{e),. . . ,Vm{s)- □ 

Clearly, literally the same arguments prove the complex analytic coun- 
terpart of this theorem. Besides, one may replace holomorphic functions by 
meromorphic (in which case the exponents in the above proof may well 
be negative). 

Remark 4. The coordinate form of Proposition 1^1 may be considered as a 
theorem on rank for meromorphic matrix functions (modulo renaming the 
parameter to a more traditional complex variable z). Any meromorphic 
germ of a n x m-matrix function X(z), z £ (C^,0), can be represented as 



X{z) = U{z)D{z)V{z), D{z) 



(3.1) 
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where the square holomorphic matrix functions U{z) and V{z) of appropri- 
ate sizes are holomorphically invertible, and D(z) is an n x m-matrix with 
only powers or zeros on the diagonal. 

Another interpretation of this result can be described as an extendability 
of holomorphic subbundles. Consider a holomorphic vector bundle over a 
one-dimension base (curve) Z. Assume that a holomorphic subbundle is 
defined over a punctured neighborhood of some point zq & Z on the base. 
If the subbundle has a "moderate" singularity, i.e., it can be spanned by 
(eventually, multivalued) sections growing no faster than polynomially at 
zq, then in fact this subbundle can be holomorphically extended at zq. 

3.2. Infinite-dimensional version. The above arguments need some 
modification to be applicable to finite-dimensional subspaces of the infinite- 
dimensional space of analytic functions of one variable. Each such subspace 
can be naturally identified with a homogeneous linear ordinary differential 
equation. 

Everywhere below d stands for the (partial) derivative d/dt; as before, 
e G (C^,0) is a one-dimensional complex parameter. 

Lemma 7. Let fi{t,e), . . . ,fn{t,£) be n functions analytic in 3 x (C"'^,0), 
which are linear independent as functions o/t G D for all e ^ 0. 

Then the functions /i, . . . , /„ satisfy a linear nth order homogeneous dif- 
ferential equation 

L^y = 0, where = ao(t, e)c?" + • • • + a„_i(t, e)c? + a„(t, e), 

with the coefficients aj analytic in D x (C^,0) and non-singular at e = 0, so 
that ao(-,0) ^ 0. 

Proof. The proof imitates that of Proposition El and goes by induction in 
n. For n = 1 the claim is obvious: the function fi(t,e) = e'^g(t,e), where 
g{-,0) ^ is analytic in B x (C^,0), satisfies the linear equation Lifi = 0, 
where Li is the linear first order differential operator, 

Li=gd-{dg), y^gy'-g'y. (3.2) 

By the choice of the coefficients of this operator depend holomorphically 
on e and the leading coefficient g{t, e) does not vanish identically at e = 0. 
The base of induction is thus established. 

Assume now that the differential operator L^-i = /i(t, e)9"~^ + ••• of 
order n — 1 with the leading coefficient /i(-, 0) 7^ 0, vanishes on the functions 
/i, . . . , /„_i, so that Ln-ifj = for all j = 1, . . . , n - 1. 

By assumption, L„_i/„ is not identically zero, otherwise /„ would be a 
linear combination of /i, . . . , fn-i for all e 7^ 0. Therefore, one can represent 
Lfi-ifn = £^g{ti£) with some finite u and 5(-,0) 7^ 0. 

Let Li be the first order operator vanishing on constructed as in (|3.2|) . 
The composition L„ = Li o L^-i vanishes on all functions /i, . . . , fn-i and 
also on /„. The coefficients of this operator are analytic functions of e and 
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the leading term is equal to g{t,e)h{t,e) which does not vanish identically 
for e = 0. □ 

For the differential operator constructed in Lemma the ratios 
Iki(') ^)ll/l|oo(') in the sense of any norm on the space of analytic func- 
tions on the disk D, are automatically bounded as e — > since the denom- 
inator does not vanish at e = 0. But since is unique modulo a common 
factor, this boundedness holds true for any analytic family of operators an- 
nulled by the functions /i, . . . , /„. 

More precisely, assume that another differential operator L* = poit, e)(?"'+ 
• • • + pn-i{t, e)d + Pn{t, e) of the same order n has polynomial coefficients 
Pj{-,e) E C[t] and defines the same set of solutions fi, ■ ■ ■ , fn for all e ^ 0. 
Denote 

b*{e) = \\p,{;e)\\c[t]>0. (3.3) 

Corollary 8. All the ratios 6^(e)/6Q(e), j = 1, . . . ,n, have finite limits as 
e ^ 0. 

Proof. Let z/j G Z, j = 0, . . . , n, be the vanishing orders of the coefficients: 
this means by definition that pj{t,e) = e'^^qj{t,e), where qj{-,e) G C[t] are 
still analytic in e, polynomial in t and qj{-,0) ^ 0. Since any two linear 
differential operators with the same null space are proportional, we have 

Pl = e-^"^o^ = ^^ j = l,...,n. (3.4) 
Po qo ao 

Since the limit values ao(-,0), qj{-,0) and qo{-,0) are all nonzero by con- 
struction, we conclude that z/j ^ fg. 

By construction, b*/bn = g'^j-^o [I'^jjjcM second ratio has a bounded 

' J' lko||c[t]' 

limit since the denominator is nonzero, and the factor e'^^~'^° is bounded. □ 

3.3. Proof of Lemma |1J Consider the polynomial family (|2.4j) satisfied 
by n entire functions analytically depending on the parameters, and as- 
sume that some of the ratios 6j(A)/6o(A) is not locally bounded. Then 
the semialgebraic set {{X,z) G i^Sn^m ^ &o(A) < zbj{X)}, a subgraph 
of the reciprocal function 6o(A)/6j(A), contains a point (A^,,0) in its clo- 
sure. By the curve selection lemma |Mil68j . one can find a real analytic 
arc 7: (M^,0) —>■ (M^" "'lA*), e 1-^ A(e), such that along this arc the ratio 
6o(A(e))/6j(A(e)) tends to zero as e ^ 0. 

But this would contradict the Corollary |S1 since the corresponding one- 
parameter family of linear ordinary differential equations L* = Lx^s) with 
polynomial coefficients, obtained by restriction of (|2.4|) on the arc 7, has an 
analytic family of entire solutions. □ 

4. Discussion and concluding remarks 

4.1. Algebra vs. analysis. Derivation of the equation ()2.4() in ^2.21 as a 

linear combination between the vectors £,0,^,1, ■■■ ,(,n over the field C(t, A) 
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is a distant analog of derivation of a linear dependence between vectors 
C, BC, B^C, S"C, where C G C" is a vector and B e Mat(n, C) a linear 
operator over the field C. However, in the latter case the coefficients of the 
linear combination are explicitly bounded. Indeed, if P{u) = + aiM"~^ + 
■ ■ ■ + ttn is the characteristic polynomial of the operator B, then P{B) = 
which translates into a linear dependence between iterates B^(^, j = 0, . . . ,n 
for any initial vector ^. The leading coefficient of this combination is 1. 
On the other hand, the magnitude of the coefficients of the characteristic 
polynomial is explicitly bounded in terms of the spectral radius of B, i.e., 
ultimately in terms of the norm \B\. Thus an analog of Lemma [5] in the 
"constant" context, is a purely algebraic fact. 

In a rather surprising way, the bound for the norms of coefficients in the 
non-constant case, is obtained using transcendental arguments, involving 
solutions of the differential equations. A challenging problem is to obtain a 
direct algebraic proof of this fact. 

4.2. Generalizations for Fuchsian systems. The condition that the sys- 
tem ()1.3p is polynomial, can be to some extent relaxed: most of the argu- 
ments would work equally well for rational systems exhibiting singularities 
(regular or irregular). 

Consider, for instance, a Fuchsian system 



with the singular points ti, . . . ,tm and the respective residue matrices Aj E 
Mat(n,C) (the point t = oo is singular if X^i^^i / 0)- The parameters 
describing this system, are n™ entries of the residue matrices, and m complex 
numbers ti,...,tm specifying the location of singular points. As before, 
denote by A G C*-" +1)™- the collection of all these data. 

Theorem 2. Any linear combination of solutions of the system 1)4. 1() satis- 
fies a linear nth order differential equation of the form 1)1.11) with polynomial 
coefficients aj G C[t], polynomially depending on the parameters A. 
// \ Aj\ ^ M , and \tj\ < M for all j = 1, . . . , m, then 



OillcM ^ llaollcM • [max(M,2)]2°*" Vj = l,...,n. (4.2) 



Sketch of the proof. The construction is completely parallel to the proof of 
Corollary El and Lemma El One has to consider the rational vector functions 
((SJ in which A{t, A) is not polynomial anymore, but rather rational matrix 
function whose entries are ratios of lattice polynomials from A]. In the 
same way as before the degrees and the norms of the numerators and the 
denominators can be explicitly controlled. The estimates will be different, 
but all difference will be finally absorbed in the 0(n^m)-term. To prove the 
Principal Lemma in this case, one should choose instead of the unit disk, 
any other disk free of singularities of the system ()4.1)1 : for each particular 
value of the parameters it is always possible. □ 




m 




(4.1) 
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Of course, Lemmas does not apply to the case when the solutions of the 
equation have singularities in the unit disk B. However, a suitable technical 
modification of this Lemma (cf. with Remark ISJ would allow to place an 
upper bound for the number of isolated zeros of solutions of (|1.H) in any 
disk free of singularities of the solutions, in terms of the size of this disk and 
the distance from it to the singular set. 

For instance, the following Corollary immediately follows from Theorem|2l 
and Remark 13 

Corollary 9. The number of zeros of any linear combination of solutions 
of the system (|4.1|) in the disk of radius r < 1 does not exceed the expression 
H1.6() provided that the concentric disk of radius 6r still contains no singular 
points of this system. □ 

Moreover, the technique developed in |RY96j and later in |NY03j allows 
to extend the upper bounds also for neighborhoods of regular singularities 
with the real spectrum and sufficiently distant from all other singularities 
of the system 1)4. Finally, one could expect on this way a significant 
improvement of the bounds given by the algorithm from |NY03j . 

4.3. Further applications. The characteristics that essentially determines 
the bounds for the number of zeros in Theorems ^ and |2J is the number n^m 
(resp., (n^ + l)m) of the "free parameters" on which depends the corre- 
sponding "universal family" of polynomial (resp., Fuchsian) systems. 

Sometimes the "universal family" appears naturally in the form depending 
on fewer number of parameters. Thus, for instance, the Picard-Fuchs system 
of differential equations for Abelian integrals has a hypergeometric form, 

{tE - B)x = Cx, B,C,E £ Mat(n, C), E = identity matrix, (4.3) 

completely determined by 2 constant matrices B, C involving altogether 
only 2n^ free parameters. The singular points of this system occur at the 
eigenvalues of the matrix B jNYOlj . 

Theorem 3. Any linear combination of solutions of the system (|4.3j) satis- 
fies a linear nth order differential equation of the form with polynomial 
coefficients aj G C[t], polynomially depending on the parameters A. 
If\B\,\C\ ^ M, then 

ll«jllcM ^ IkollcM • [max(M, 2)]^°*" Vj = l,...,n. (4.4) 

Proof. The proof is completely similar to the proof of Theorems Hand 121 □ 

Note that this is a better bound than the one obtained by the straight- 
forward reduction of the system (|4.3j) to the rational system 

A(t) 

X = ^^^^^^ -B) '^' eMat(n,C[t]), degA^n, 

involving n'^ + n free parameters. Of course, an analog of Corollaryl^is valid 
in the hypergeometric case as well. 
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